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Abstract 

We introduce the first analytical model of asymmetric community dynamics 
to yield Hubbell's neutral theory in the limit of functional equivalence among 
all species. Our focus centers on an asymmetric extension of Hubbell's lo- 
cal community dynamics, while an analogous extension of Hubbell's meta- 
community dynamics is deferred to an appendix. We find that mass-effects 
may facilitate coexistence in asymmetric local communities and generate uni- 
modal species abundance distributions indistinguishable from those of sym- 
metric communities. Multiple modes, however, only arise from asymmetric 
processes and provide a strong indication of non-neutral dynamics. Although 
the exact stationary distributions of fully asymmetric communities must be 
calculated numerically, we derive approximate sampling distributions for the 
general case and for nearly neutral communities where symmetry is broken by 
a single species distinct from all others in ecological fitness and dispersal abil- 
ity. In the latter case, our approximate distributions are fully normalized, 
and novel asymptotic expansions of the required hypergeometric functions 
are provided to make evaluations tractable for large communities. Employ- 
ing these results in a Bayesian analysis may provide a novel statistical test to 
assess the consistency of species abundance data with the neutral hypothesis. 
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mass-effects 



1. Introduction 



The ecological symmetry of trophically similar species forms the central 
assumption in Hubbell's unified neutral theory of biodiversity and biogeog- 

In the absence of stable coexistence mechanisms, 

a stochastic pro- 



raphy (Hubbell 



200ip . 

local communities evolve under zero-sum ecological drift 
cess of density-dependent birth, death, and migration that maintains a fixed 
community size (Hubbell, 2001). Despite a homogeneous environment, mi- 
gration inhibits the dominance of any single species and fosters high levels of 
diversity. The symmetry assumption has allowed for considerable analytical 
developments that draw on the mathematics of neutral population genet- 



ics (] 


^isher 


1930 Wright 


1931) to derive exact predictions for emergent, 


macro-ecological patterns ( 


Chave 2004| E 


Etienne and Alonso} 


2007| 


McKane 


et al. 


2000 


Vallade and Houchmandzadeli 


. 2003 


Volkov et al. 


. 2003 


Etienne 


and Olffl 20041 McKane et al. 2004[ Pigolotti et al., 2004[ Hel |2005; 


Volkov 


et al. 


2005; 


Hu et al. 2007 


[ Volkov et al. , 


2007 


Babak and He] [2008 


2009). 



Among the most significant contributions are calculations of multivariate 
sampling distributions that relate local abundances to those in the regional 



enne 



metacommunity (Alonso and McKane, 2004; Etienne and Alonso, 2005; Eti 



2005, 2007). Hubbell first emphasized the ultility of sampling theories 
for testing neutral theory against observed species abundance distributions 
(SADs) (Hubbell, 2001). Since then, Etienne and Olff have incorporated 



sampling distributions as conditional likelihoods in Bayesian analyses (Eti- 



enne and Olff, 2004, 2005 Etienne, 2007, 2009). Recent work has shown that 



the sampling distributions of neutral theory remain invariant when the re- 



striction of zero-sum dynamics is lifted (Etienne et al. 2007 Haegeman and 



Etienne, 2008 Conlisk et al. 2010) and when the assumption of strict sym- 



metry is relaxed to a requirement of ecological equivalence (Etienne et al. 



2007 


Haegeman and Etienne 


et al. 


2009 


)• 



The success of neutral theory in fitting empirical patterns of biodiver- 



sity (Hubbell, 2001 Volkov et al., 2003, 2005 He, 2005 Chave et al., 2006) 



has generated a heated debate among ecologists, as there is strong evidence 



for species asymmetry in the field ( 


Harpei 


1977 


Goldberg and Barton 


1992 


Chase and Leibold 


2003 


Wootton 


2009 


1/evine and HilleRisLambers 


2009) 



Echoing previous work on the difficulty of resolving competitive dynamics 
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from the essentially static observations of co-occurence data (Hastings, 1987), 
recent studies indicate that interspecific tradeoffs may generate unimodal 



SADs indistinguishable from the expectations of neutral theory (Chave et al. 



2002 



2005 



Mouquet and Loreau 2003| Chase[ |2005| |He[ |2005| |Purves and Pacala 



Walker, 2007; Doncaster, 2009). These results underscore the compat- 



ibility of asymmetries and coexistence. The pioneering work of Hutchinson 



(1951), has inspired a large literature on asymmetries in dispersal ability 
that permit the coexistence of "fugitive species" with dominant competitors. 
In particular, Shmida and Wilson (1985) extended the work of Brown and 



Kodric-Brown (1977) by introducing the paradigm of "mass-effects", where 
immigration facilitates the establishment of species in sites where they would 
otherwise be competitively excluded. Numerous attempts have been made 
to reconcile such deterministic approaches to the coexistence of asymmetric 
species with the stochastic model of ecological drift in symmetric neutral 



theory ( 


Zhang and Lin 


19971 Tilman, 


2004; 


Chase 


2005 


Alonso et al. 


Gravel et al. 


|2006[ |Pueyo et al-l |2007| [Walker 




>007| 


Alonso et al. 


Ernest et al. 
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. Many of these attempts 



2006 



2008 



on insights from the concluding chapter of Hubbell's book (Hubbell, 2001). 
Nevertheless, the need remains for a fully asymmetric, analytical, sam- 



pling theory that contains Hubbell's model as a limiting case (Alonso et al. 



2006). In this article, we develop such a theory for local, dispersal-limited 



communities in the main text and defer an analogous treatment of meta- 
communities to Appendix A. Hubbell's assumption of zero-sum dynamics is 
preserved, but the requirement of per capita ecological equivalence among 
all species is eliminated. Asymmetries are introduced by allowing for the 
variations in ecological fitness and dispersal ability that may arise in a het- 
erogeneous environment (Leibold et al. , 2004 Holyoak et al. 2005 ). Our work 
expands on the numerical simulations of Zhou and Zhang (2008), where vari- 
ations in ecological fitness alone were considered. Coexistence emerges from 
mass-effects as well as ecological equivalence, and both mechanisms generate 
unimodal SADs that may be indistinguishable. For local communities and 
metacommunities, we derive approximate sampling distributions for both the 
general case and the nearly neutral case, where symmetry is broken by a sin- 
gle species unique in ecological function. These approximations yield the 
sampling distributions of Hubbell's neutral model in the limit of functional 
equivalence among all species. 
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2. A general sampling theory for local communities 



For a local community of Jl individuals and S possible species, we model 
community dynamics as a stochastic process, N(r), over the labelled commu- 
nity abundance vectors n = (n 1; . . . , rig). Consistent with zero-sum dynam- 
ics, we require all accessible states to contain Jl total individuals: 52 f=1 rii = 
Jl and < < Jl • The number of accessible states is A = J2ni=o ' ' ' Sns=o 
5{J L -ni n s ). 

Allowed transitions first remove an individual from species i and then 
add an individual to species j. Removals are due to death or emigration 
and occur with the density-dependent probability m/ Jl- Additions are due 
either to an immigration event, with probability mj, or a birth event, with 
probability 1 — rrij. We will refer to the mj as dispersal abilities. If im- 
migration occurs, we assume that metacommunity relative abundance, Xj, 
determines the proportional representation of species j in the propagule rain 
and that the probability of establishment is weighted by ecological fitness, 
Wj, where high values correspond to a local competitive advantage or a supe- 
rior adaptation to the local environment. Therefore, species j recruits with 
probability 

(1) 



w jXj 



2^k=l W k x k 

where Xj G (0,1), Wj G (0, oo), and ^2 k=1 Xk = 1. If immigration does not 
occur, we assume that local relative abundance, uj/Jl, governs propagule 
rain composition such that species j recruits with probability 



,s — ■ ( 2 ) 



J2 k=1 w k n k - Wi 



In numerical simulations of an asymmetric community, Zhou and Zhang 



(2008) employed a similar probability for recruitment in the absence of immi- 
gration. Here, a factor of Wi is subtracted in the denominator because species 
i loses an individual prior to the birth event for species j. An analogous sub- 
traction is absent from Eq. [T] because we assume an infinite metacommunity 
where the Xj are invariant to fluctuations in the finite, local community pop- 
ulations. 
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In sum, the nonzero transition probabilities are stationary and given by 




Pr{N(r + At) = n - e { + eAN(T) = n} 
jH = &$o At 

- mj) WjHj + m s s pj ) . (3) 

2^k=l W k n k - Wi 2^k=l W k x k / 

where e« is an S'-dimensional unit vector along the ith-direction, the Wk 
must be sufficiently large such that X/fc=i w k n k — Wi > 0, and the time, r, is 
dimensionless with a scale set by the overall transition rate. The probability 
of state occupancy, Pft, evolves according to the master equation 

dP- S S 

-~T~- — {Tijn+St-ejPn+et-ej — TjirfPfi) Qij, (4) 

i=l j=l,j^i 

where 

Qij = e(J L -(^ + l))6(n J -l), (5) 

and we define the step-function 0(x) to be zero for x < and one otherwise. 
Eq. [4] can be recast in terms of a transition probability matrix W 

rip A 

rE* ( 6 ) 

6=1 

where a, b e (1, . . . , A) enumerate accessible states with components (ai, . . . , 
ag), (61,..., 65). The left eigenvector of W with zero eigenvalue yields 
the stationary distribution for community composition, P* = lim^oo P a {^)- 
Marginal distributions yield the equilibrium abundance probabilities for each 
species i 

A 

0=1 

From here, we calculate the stationary SAD by following the general treat- 
ment of asymmetric communities in Alonso et al. (2008) 

s 

s: = J2 p n ] *- ( 8 ) 

i=l 
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The expected species richness is 



S* = J2 S n <S - 



(9) 



n=l 



Given that the local community, with abundances nj, is defined as a sample 
of the metacommunity, with relative abundances Xi, we have established the 
framework for a general sampling theory of local communities. 

This sampling theory incorporates aspects of the mass-effects paradigm 



(Brown and Kodric-Brown 


1977 


Shmida and Wilson 


1985; 


Holt 


1993 


Lei- 


bold et al. 


2004; 


Holyoak et al. 


2005 


). Local asymmetries in ecological fit- 



ness imply environmental heterogeneity across the metacommunity such that 
competitive ability peaks in the local communities where biotic and abiotic 



Tilman 


1982 


Leibold , 


1998 



Chase and Leibold, 2003). Where species experience a competitive disad- 
vantage, the mass-effects of immigration allow for persistence. Indeed, the 
master equation given by Eq. |4| when applied to open communities where 
rrij > for all j, admits no absorbing states and ensures that every species 
has a nonzero probability of being present under equilibrium conditions. By 
contrast, when Eq. [4] is applied to closed communities where rrij = for all 
j, the eventual dominance of a single species is guaranteed. 

Mass-effects allow for a soft breaking of the symmetry of neutral theory 
and provide a mechanism for multi-species coexistence. In Fig. [TJ we present 
numerical results for the marginal equilibrium distributions of an asymmetric 
local community subsidized by a potentially neutral metacommunity, where 
the five species share a common relative abundance, Xj=0.2. Although a sin- 
gle species may dominate due to a locally superior competitive ability (see 
Fig. [jji), multi-species coexistence may arise, despite significant competitive 
asymmetries, due to high levels of immigration that tend to align local rel- 
ative abundances with those in the metacommunity (see Fig. [lj). Despite 
the underlying asymmetric process, coexistence via mass-effects generates 
unimodal SADs that, given sampling errors in field data, may be indistin- 
guishable from SADs due to neutral dynamics, as shown in Fig. [l]i. This 
reinforces previous conclusions that the static, aggregate data in unimodal 
SADs cannot resolve the individual-level rules of engagement governing the 



origin and maintenance of biodiversity ( Chave et al. , 2002 



Mouquet and 



Loreauj [2003j |Purves and Pacala[ |2005| |He] [2005] [Chasej [2005 



Walker[ |2007] 



Doncaster, 2009). However, SADs with multiple modes are not uncommon in 
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Figure 1: Equilibrium abundance probabilities and corresponding SADs for an asymmetric 
community of Jl = 75 individuals and S = 5 species. Non-neutral dynamics generate local 
deviations from the relative metacommunity abundance, Xi = 0.20, common to all species, 
(a) Dominance by the species with highest ecological fitness given competitive asymmetries 
in a community of uniform dispersal abilities (the m, = 0.10 for all i). (b) The resulting 
bimodal SAD provides a strong indicator of non-neutral dynamics, (c) Coexistence arising 
from mass-effects (the nii — 0.90 for all i). (d) The resulting unimodal SAD closely 
resembles the SAD for a symmetric community of Jl = 75 individuals and 5 = 5 species 
where the m, — 0.35 and the x,- = 0.20. 
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nature (Dornelas and Connolly, 2008 Gray et al. , 2005) and provide a strong 
indicator of non-neutral dynamics (Alonso et al. 2008). Fig. [TJd presents a 
bimodal SAD for an asymmetric local community with low levels of immi- 
gration. 

Each plot in Fig. [T] displays results for a relatively small community of 
J L = 75 individuals and 5 = 5 possible species. Sparse matrix methods 
were used to calculate the left eigenvector with zero eigenvalue for transition 
matrices of rank ~ 1.5M. Obtaining stationary distributions for larger, more 
realistic communities poses a formidable numerical challenge. This motivates 
a search for analytically tractable approximations to sampling distributions 
of the general theory. 



3. An approximation to the sampling distribution 



The distribution Pi is stationary under Eq. [4] if it satisfies the condition 
of detailed balance 

Tijn+ei-SjPfi+ei-ej = ^jin^ni (10) 

for all i and j such that i ^ j and 6^ 7^ 0. For general (g) large- Ji 
communities where S, Wk « Xw=i w i n i l° r an we will show that detailed 
balance is approximately satisfied by 

P? = Z ? L J \ ) f[ < fe (l - m k r (M nk , (11) 



where 



j (Ef=i«w - w fe )/(j L - 1) 
Li=i w i x i 



(12) 



where (a) n = T(a + n)/T(a) is the Pochhammer symbol, Z~ l is a normal- 
ization constant, and Ik = rnk(JL — 1)/(1 — m k) is a generalization of the 
"fundamental dispersal number" (Etienne and Alonso, 2005). From the def- 
inition of Tijfi in Eq. [3j we have 

Tjjn+gj-gj rij + l wj l- rrij n j -l + 4> jH Xj 

rp — ; _ / ...... \ 1 v 1 ^/ 

i«™ n i Wi 1 m * rii + 4>iHXi 1 1 + =^ 



J2i=i wini-Wi 
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and assuming the form of PS* in Eq. 1 1 , we find 



P|* Hi + 1 Wj 1 — -p-r (^fcn^fc), 

m . i/i. 1 <m . X -L 



P ft+Si-S j n i w * 1 m * fe= i ( 0^x fc ( 1 + — 



(14) 

Now, for large- Jl communities where Wk « z2i=i w i n i f° r au ^? the ratio 
e i:; fc = (tUi — Wj)/(^2 l=l Wini — Wk) is a small number. Given (a(l + e)) n ~ 
(a) n + (9(e), we expand the right-hand-side of Eq. 14 to obtain 

pg* _, _ s 



PS * T ■— 



+ J20(e ljk ) } (15) 



which validates our assertion that Eq. fTTl is an approximate sampling dis- 
tribution of the general theory when S « Y2i=i w i n i- For communities of 
species that are symmetric (s) in ecological fitness but asymmetric in disper- 
sal ability, Eq. 11 reduces to an exact sampling distribution 

S 

P « = Z "{nJ L .^n)^- mk)nk{IkXk) -^ (16) 

that satisfies detailed balance without approximation. Analogous distribu- 
tions for general and fitness-symmetric metacommunities are provided in 
Appendix A. However, in all of these results, the normalization constants 
must be calculated numerically. This limits the utility of our sampling dis- 
tributions in statistical analyses. Can we find a non-neutral scenario that 
admits an approximate sampling distribution with an analytical expression 
for the normalization? 



4. Sampling nearly neutral communities 

As the species abundance vector evolves under Eq. [4j consider the dynam- 
ics of marginal abundance probabilities for a single focal species that deviates 
in ecological function from the surrounding, otherwise symmetric, commu- 
nity. In particular, let the first element of N(r) be the marginal process, 
N(t), over states n € (0, . . . , Jl), for the abundance of an asymmetric focal 
species with dispersal ability m, ecological fitness w, and relative metacom- 
munity abundance x. If all other species share a common dispersal ability 



9 



m and ecological fitness w Q , then the focal species gains an individual with 
probability 

* - E T «(« 

i=2 

Jl — n ( x urn to 



1 — TO) 77; r + TO 



Jl \ urn + w Q ( Jl — n — 1) to + - 1) 

and loses an individual with probability 

i=2 

n ( w (J L -n) w a (l-x) 
I - m )— — — + m 



(17) 



Jl V w{n — 1) + w (Jl — n) wx + w (l — x) 



(18) 



where we have used Ylt=i x k — 1- These marginal transition probabilities do 
not depend separately on w and w a , but only on their ratio. Without loss of 
generality, we redefine w = w/w Q to be the focal species' local advantage in 
ecological fitness. Eqs. 17 and 18, which are independent of the abundances 
(n 2 , . . . ,ns), suggest a univariate birth-death process for the marginal dy- 
namics of the asymmetric species governed by the master equation 

dP 

— ^ = g n -Mn - l)P n _! + r n+1 9( J L - (n + l))P n+1 

GST 

-{g n Q{JL - (n + 1)) + r n Q(n - l))P n , (19) 

and we formally derive this result from Eq. [4] in Appendix B. Given the 
well-known stationary distribution of Eq. [19] 



n-l 

9 



P * n (20) 



we find an exact result for the stationary abundance probabilities of the focal 
species in a nearly neutral (nn) community 

pnn* _ z (Jl\ B(\ + n,tj-n) 

Pn ~ Z \n) V B(A,0 ' (21) 

10 



where B(a, b) = T(a)T(b)/T(a + b) is the beta-function 

Z- 1 = 2 F 1 (-J L ,X;1-^, V ), 



(22) 



and 



7/ 

A 



w- 



1 — m + x(w — 1) 
1 — m u; + x(w — 1) ' 
(Jl — l)mx 
1 — m + x(w — 1) ' 

(Jl — 1)(1 — xwm + x(u? 



1 + 



1)) 



(23) 



1 — wm + x(w — 1) 

For the asymmetric focal species, this is an exact result of the general model, 
Eq. |4j that holds for nearly neutral local communites with any number of 
additional species. Eq. 21_ may be classified broadly as a generalized hyperge- 
ometric distribution or more specifically as an exponentially weighted Polya 
distribution (Kemp, 1968 Johnson et al. , 1992). 



In the absence of dispersal limitation, Eq. 21 becomes 

Jl 



lim PI 

m,m — >1 



1 + x(w — 1) 



Jl 

n 



[vox 



) n (l-x) jL - n , 



(24) 



where the identity B(a, b)B(a + b, 1 — b) = ir/(a sin(7r6)) has been used. This 
is a weighted binomial distribution with expected abundance wxJl/(1 + 
x(w — 1)) and variance wx(l — x)Jl/(1 + x(w — l)) 2 . In the neutral, or 
symmetric, limit where w — 1, Eq. 24 reduces to a binomial sampling of the 
metacommunity, sensu Etienne and Alonso (2005). 

In the presence of dispersal limitation, we evaluate E^nP" 11 * to obtain 
the expected abundance 



E[N*} 



9 1 7 

JlXv 2^1 (1 



J L ,l + A;2-£;?7) 



where iV* 
given by 



= lim T 



Var[A^* 



e-i 2F!(-j l ,x-a-^v) ' (25) 

. The variance of the stationary distribution is 



d d 
or] or] 
E[N* 2 ] - E[iV*] 2 , 



(26) 
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and we evaluate S^ 1 n 2 P° n * to obtain 



E[N 



*2l 



J L X V 3 F 2 (1 - J L , 1 + A, 2; 2 - f, 1; rj) 



^-1 



2 F 1 (-J i ,A;l-£; 77) 



(27) 



In Eqs. 25 and 26 the normalization of Eq. 22 generates central moments for 
the abundance distribution and plays a role analogous to the grand partition 
function of statistical physics. Recent studies have demonstrated the utility 



of partition functions in extensions of Hubbell's neutral theory (O'Dwyer 



et al. 2009 O'Dwyer and Green, 2010) 



For large- Jl communities, evaluation of the hypergeometric functions in 
Eqs. [211 [25l and [27] is computationally expensive. To remove this barrier, 



one of us (N.M.T.) has derived novel asymptotic expansions (see Appendix 
C). We use these expansions to plot the stationary abundance probabili- 
ties for Jl = 1M. In Fig. [2^,, small local advantages in ecological fitness 
generate substantial increases in expected abundance over the neutral pre- 
diction. Hubbell found evidence for these discrepancies in Manu forest data 
and referred to them as "ecological dominance deviations" ( jHubbell 2001). 
Hubbell also anticipated that dispersal effects would mitigate advantages in 
ecological fitness (Hubbell, 2001). The right panel of Fig. [2] demonstrates, 



once again, that enhanced mass-effects due to increased dispersal ability may 
inhibit the dominance of a locally superior competitor by compelling relative 
local abundance to align with relative metacommunity abundance. 

An approximation to the multivariate sampling distribution of nearly 
neutral local communities is constructed in Appendix B 



z| Jl \ n B(A + n,£-n 
K n,n 2 , . . . ,n s ) 



B(A,0 ((1-x) 



'on) J^—n 



IK 

i=2 



where 



1 + n{w - l)/( J L - 1) 



'on ± o 



(28) 

l + x(w-l) ' ^ 

A related approximation for the sampling distribution of nearly neutral meta- 
communities is derived in Appendix A. In the absence of dispersal limitation, 
Eq. [28] becomes 



lim PI 

m,m — >1 



1 



Jl 



1 + x(w - 1) J \n,n 2 ,---,ns 



Ji 



n 

i=2 



:<rx) n I I x™\ 



(30) 
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Figure 2: Novel asymptotic expansions of hypergeometric functions have been used to 
plot marginal equilibrium abundance probabilities for the asymmetric focal species, with 
relative metacommunity abundance x — 0.01, in a nearly neutral local community of 
J/,=1M individuals, (a) Dominance with rising advantage in ecological fitness, as indicated 
for each curve. Here, all species are symmetric in their dispersal ability (m = m = 0.10). 
(b) Dispersal mitigates the advantage in ecological fitness (w = 1.08) of the asymmetric 
focal species. All other species share a common dispersal ability of m — 0.10. 
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where we have used (a) n ~ a n + 0(a n 1 ) for large a. Finally, in the symmetric 



limit, Eq. [30] reduces to a simple multinomial sampling of the metacommu- 
nity, as expected. 

To illustrate the impacts of an asymmetric species on the diversity of 
an otherwise symmetric local community, Fig. [3] plots Shannon's Index of 
diversity 

// = _MQ log Em_f B log B, (31) 

1=1 

for various values of the ecological fitness advantage, w, and dispersal ability, 
m, in a nearly neutral community of S = 5 species and Jl = 75 individuals. 
All five species share a common relative metacommunity abundance, x = 



Xi = 0.2, so given the exact result for E[iV*] in Eq. 25, we know immediately 
that E[iV*] = (Jl — E[N*])/(S — 1) for the remaining symmetric species. 
Note that H is maximized where all abundances are equivalent, such that 
E[N*]/J L = E[N*]/J L = Xi. As can be seen from the next section, this 
relation holds in the neutral limit where w = 1 and m = m = 0.1, but 
small asymmetries in dispersal ability have a negligible impact on diversity 
when all species are symmetric in ecological fitness. Therefore, each curve 
in Fig. [3] peaks near w — 1 at approximately the same value of H. Away 
from w = 1, the declines in diversity are regulated by mass-effects, with more 
gradual declines at higher values of m. 

5. Recovering the sampling distribution of neutral theory 

In a perfectly symmetric local community, the stochastic dynamics for 
each species differ solely due to variations in relative metacommunity abun- 
dances, the Xi. In particular, if rrij = m and Wj = 1 for all j in Eq. |3j we 
recover the multivariate transition probabilities for a neutral sampling the- 



ory of local communities, as suggested on p. 287 of Hubbell's book (Hubbell 



2001). Similarly, in the symmetric limit of Eq. 19 where m = m and w = 1 



we recover the marginal dynamics for neutral (n) local communities with 



stationary distribution (McKane et al. 2000) 



J L \B(Ix + n,J L + 1(1 - x) - n) 

5 



n J B(Ix,I(l -x)) 

where / = m(J l — 1) / (1 — m) . This result follows from the symmetric limit of 



Eq. 21 after applying the identity r(a)r(l — a) = it/ sin(7ra). The expected 
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Figure 3: Plots of the information-theoretic diversity metric, H , for a nearly neutral local 
community of Jl = 75 individuals and 5 = 5 species over various values of the ecological 
fitness advantage, w, and dispersal ability, m. The symmetric species share a common 
dispersal ability of m — 0.1, and all species share a common relative metacommunity 
abundance such that x = xi = 0.2. Diversity peaks where expected local abundances are 
equivalent, and this occurs in the symmetric limit, given by w = 1 and m = m . The 
asymmetries in dispersal ability shown here have a negligible impact on diversity when all 
species are symmetric in ecological fitness, so each curve in Fig. [3] peaks at approximately 
the same value of H near w = 1. Away from the peak, declines in diversity are regulated 
by mass-effects, with more gradual declines at higher values of m. 
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abundance and variance are obtained from the symmetric limits of Eqs. 25 
and [26} respectively, after applying the identities in Eqs. C. 1.0.1 and C.2.2.2 



E[iV*] 
Var[iV*] 



xJ L , 

x(l — x)Ji 



Jl + I 
1 + / ' 



(33) 
(34) 



Finally, the symmetric limits of Eqs. 11, 16, and 28 all yield the stationary 



sampling distribution for a neutral local community (Etienne and Alonso 



2005 Etienne et al., 2007) 



J I 



m, . . .,nsj (I) 



Ml 



(35) 



In the special case of complete neutrality, Eq. 35 is an exact result of the 
general model, Eq. [4j This sampling distribution continues to hold when 



the assumptions of zero-sum dynamics and stationarity are relaxed (Etienne 



et al. 2007 Haegeman and Etienne 2008). 



6. Discussion 



We have developed a general sampling theory that extends Hubbell's 
neutral theory of local communities and metacommunities to include asym- 
metries in ecological fitness and dispersal ability. We anticipate that a pa- 
rameterization of additional biological complexity, such as asymmetries in 
survivorship probabilities or differences between the establishment probabil- 
ities of local reproduction and immigration, may be incorporated without 
significant changes to the structure of our analytical results. Although the 
machinery is significantly more complicated for asymmetric theories than 
their symmetric counterparts, some analytical calculations remain tractable. 
We find approximate sampling distributions for general and nearly neutral 
communities that yield Hubbell's theory in the symmetric limit. Our fully 
normalized approximation in the nearly neutral case may provide a valuable 
statistical tool for determining the degree to which an observed SAD is con- 
sistent with the assumption of complete neutrality. To facilitate a Bayesian 
analysis, we have enabled rapid computation of the required hypergeometric 
functions by deriving previously unknown asymptotic expansions. 
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Appendix A. Sampling asymmetric metacommunities 



The analytical insights of Etienne et al. (2007) suggest a clear prescription 



for translating local community dynamics into metacommunities dynamics 
in the context of Hubbell's unified neutral theory of biodiversity and bio- 



geography (Hubbell, 2001): replace probabilities of immigration, rrij, with 
probabilities of speciation, Vj\ assume Xj ~ I/St + 0(1/ S T ) for all j, where 
St is the total number of species that could possibly appear through speci- 
ation events; and consider asymptotics as St becomes large. 

Following this recipe, we translate the transition probabilities for asym- 
metric local communities, Eq. [3j into the transition probabilities for asym- 
metric metacommunities (M) 



M 
ijn 




1 - V. 



WjUj 



3 ) st—\St 



+ V 



Wi 




EST 
k=l W k 



(A.l) 



where Jm is the number of individuals in the metacommunity, Wj/ J2k=i w k 
is the probability that an individual of species j establishes following a spe- 
ciation event, and 



ej f = 0(J M -(n i + l))6(n i -l). 
Metacommunity dynamics are governed by the master equation 



(A.2) 



dr 



St Sjp 

E E 

i=l j=l,j^i 



pM rpM pM \ r\M 

n+ei—Sj ii+ei—ej jin n J ij 



iAf r>M 



M 



(A.3) 



If Uj > for all j, there are no absorbing states, so for large- St, there is 
a nonzero probability that any given species j exists. Analogous develops 
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to those in Section [3] show that detailed balance in the general theory is 
approximated by 

p * M ' = z ^L J . M ns )n ro '' (i -*)"*(*)»■ < a - 4 > 



where 

iM fl (E£'l^-W fc )/(JM-1) s 
0fcn = ^ ^SV ' l A ' 5 J 

and 6k = Vk(Ju — 1)1(1 — ^fc) is the generalization of Hubbell's "fundamental 
biodiversity number" (Hubbell |2001 ). The fitness-symmetric (s) distribution 



P^ = Z~d Ju ) fl(l - v k r (0 k /S T ) nk , (A.6) 

satisfies detailed balance up to 0(1/ St)- 

For the special case of nearly neutral metacommunities, we translate the 
marginal dynamics for an asymmetric species in an otherwise symmetric 
local community into the marginal dynamics for an asymmetric species in an 
otherwise symmetric metacommunity. The transition probabilities are 

m Jm-u ( wn w ( 1 



m n ( n \ Jm-u ( w\ ( 1 

J M V j m + n(w - 1) - w V s t J \S T 



w 
St 


+ o 


1 - 


w \ 


St) 



(A.7) 



where the asymmetric focal species has speciation probability v and enjoys an 
ecological fitness advantage, w, over all other species, which share a common 
probability of speciation v . If n = is an accessible state, then as St 
becomes large and w remains finite, the equilibrium probability of observing 
the asymmetric species approaches zero. However, if we assume that the 
asymmetric species is identified and known to exist at nonzero abundance 
levels, the stationary distribution is 

pnn,M* y ( J m\ n B(A M + n,£ M ~ ft) / a ^ 

Pn = Zm { n P M ' (A - 8) 
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with 
and 



Z K l = 2-Fi(- Jm, Am; 1 - £m; Vm) - 1, (A.9) 



(J M - 1)(J M + 9o - 1) n / J_ 

%/ " W (Jm-^o(^-1)-1)(J A / + ^-1) Ut 



9 c/ 1 



C 1 ~ \ C2 



= 1)( Jm + ,„ 1) + /IN 

•/Af - o [w - 1) - 1 \S T J 

where 6 = (Jm — 1)^/(1 — z/ ) and O = (Jm — l)v /(l — z/ Q ) are Hubbell's 
"fundamental biodiversity numbers" for the asymmetric species and all other 
species, respectively. 

An approximate multivariate stationary distribution is obtained in an 



identical manner to the derivation of Eq. 28 



D nn,M* rv I -Jm \ n B(A A f + 71, £m ~ 7i) 



K n,n 2 , . . . ,n ST J B(A M ,6v/ 
1 

x 



((1 - l/S T )(j) n>M )j u . 



H^m/StU, (A.ll) 



where 



1 + 1^) +0 (^- (A12) 



We now propose a modest extension to the prescription in Etienne et al.| 



(2007) for converting multivariate distributions over labelled abundance vec- 
tors to distributions over unlabelled abundance vectors. Because the asym- 
metric focal species has been identified and is known to exist with abundance 
n > 0, this species must be labelled, while all other species are equivalent 
and may be unlabelled. Therefore, we aim to transform Eq. A.ll into a mul- 
tivariate distribution over the "mostly unlabelled" states n = (n, n 2 , . . . , ns), 
where 5 is the number of species observed in a sample and each (n 2 , ■ ■ ■ , ns) 
is an integer partition of Jm — n. (To provide an example, if Jm = 3, four 
distinct states are accessible: (3) with 5=1, (2,1) with 5 = 2, (1,2) with 
5 = 2, and (1, 1, 1) with 5 = 3.) The conversion is given by 



(S T - 1) ■_ pim,M* /a -|q\ 

lli=0 
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where is the number of elements in (n 2 ,. . . ,hs) equal to i. Note that 
$o = St — 1 — (S — 1). Taking the leading behavior for large- St, we obtain a 
modification of the Ewens (1972) sampling distribution appropriate to nearly 



neutral metacommunities 

{s T -i)\ 



lim P. 

St— too 



nn.M* 



lim 

St ^ (s T -s)\ n nf =2 fii n£r n ^ 

V n M (d/S T )n^M)-n 



X ■ 



X 



n - l)!( 2 Fi(- J M , 1 - 6*; - 1) 



((1 - l/5 , T )0 niM )j A/ -n 



n 

i=2 



(0n,M/ •S't)^ 



n 5 - 



1 ! 



(53 



5t->oo — if))! \^SV 



J 



M ■ 



X 



2 F 1 (-J M ,e/S T ; 1 - Cm; »7Af) - 



lS-1 



M —n , 



n(J M ~n) 



n) 



'n,M ) Jm—n 

b s ~ l 



M )—n 



nf= 2 ^n/=r" 



(A.14) 



where (a) = 1 allows us to take a product over the observed species, S, 
rather than the total number of possible species, St, in the first expression; 
(z) n /(n — 1)! ~ z + 0(z 2 ) as z approaches for n > has been used to obtain 
the second expression; l'Hopital's rule along with lim^o ^2-^1(0, b; c; z)/db = 
az 3 F 2 (a + 1, 1, 1; c + 1, 2; z)/c has been used to obtain the third expression; 
and 

Zm = ^73^2(1 - Jm, 1, 1; 2 - Cm, 2; r/ M ), (A.15) 

with asymptotics of the hypergeometric function provided in §Appendix C.3| 
and §Appendix C.4 In the neutral limit, we obtain a modification to the 
Ewens sampling distribution for the scenario where a single species is labelled 
and guaranteed to exist 



lim P 

St— too 



n,M* 



e 



J 



M 



J, 



M ■ 



6 s 



e+i e+j 



M 



(A. 16) 
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Converting this result to a distribution over the "fully unlabelled" states 
n = ( fix , . . . , h s ) , we multiply by 



Jm 



e + i e + j j 



M 



(A. 17) 



and recover the Ewens sampling distribution ( |Ewens 1972), which is also the 
sampling distribution for Hubbell's metacommunity theory (Hubbell, 2001). 



Appendix B. Marginal dynamics for the local community 

We first demonstrate that the marginal dynamics of the asymmetric 
species in Eq. [19] can be derived from the multivariate dynamics of Eq. |4j 
Let 

Y,6= J2...J25(J L -n-n 2 n s ), (B.l) 

712=0 ns=0 

so that the marginal distribution for the asymmetric species is given by 



J L -n 
P n = J2 5P n- 



(B.2) 



Applying Eq. | B.1| to both sides of Eq. [4], we obtain 

dPn 



dr 



8=1 j = l,j^l 

J L -n S S 

E 6 H E (T ijH+ ^ gj P n+ ^ gj - TjiftP.fi) 6 4 

i=2 j=2,j^i 
J L -n S 

+ ^ $ (^lj'rt+ei — e"j -f n+ei — e*j ~~ TjlftPft) @lj 

i=2 

i=2 



(B.3) 
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By inspection, the first term is identically zero and the remaining terms 
generate the right-hand side of Eq. 19, namely 

r n +MJ L - (n + l))P n+l , 



J L -n S 

2j 8 T jlfi P H Q 1:j 

J L -n S 

8 Tjifj + g i -g 1 Pfi + g i -g 1 Qii = g n -iQ(n — l)P n _i 



^6(J L -(n + l))P n , 



i=2 



Jz,-n 5 



r n e(n - l)P n . 



(B.4) 



i=2 



To provide an illustration, let Jl = 4, 5 = 3, and n = 1. Then, 

J L ~n S 

— (P(2,0,2)(Pl2(2,0,2) + Pl3(2,0,2)) + P(2,l,l) (^12(2,1,1) + ^13(2,1,1)) 

+P(2,2,0)(Pl2(2,2,0) + ^13(2,2,0))) ~ (n + 1)) 

= r n+1 0(J L - (n + 1))(P(2A2) + P(2,l,l) + P(2,2,0)) 

= r n+1 6(J L -(n + l))P„, +1 , (B.5) 
where we have used the definitions of 0jj from Eq. [5j r„ from Eq. 18 , and 



P n from Eq. B.2 



We construct an approximation to the multivariate sampling distribution 
of a nearly neutral community, P? n *, by following the subsample approach 
of Etienne and Alonso (2005) and Etienne et al. (2007) that centers on the 

(B.6) 



identity. 



5-1 

pnn* pnn* "1 T pnn* 

n n J_ J_ nf\n,n2,...,n,f_i' 

/=2 



Assuming that P n , P n2 \ n , . . . , Pn f _ 1 \n,n 2 ,...,n f _ 2 are nonzero and stationary, we 
argue that conditional marginal dynamics for P n/ |n,n 2 ,...,n / _ 1 are approximated 
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by the master equation 

dP nf \n,n 2 ,...,n f _ 1 £\( _l\p 

+I|n,n2,...,n^_ 1 

- (<& v |fi ) n 3 ,...,» / _i.6(«JL ~ { n f + 1)) 

"I"' r n^|n,n2,...,n^_ 1 ©(^Z 1)) Pnf\n,n2,...,rif-n (-^-7) 

where 

„K» .-t^ . «, 

i=/+l 

= ^-% / \ 

J L [} m °)j L + n(w-l)-l +m °l + x(w-l))> 

^ ^t^ m „ „ s) 

t=/+l 



Jl V 1 J ^ + n(w - 1) - 1 + m °l + s( w - 1) )> 



(B.8) 



and 



/ 

k=2 
f 

Xf = x + ^Xfc. (B.9) 

k=2 

Eq. B.7 can be derived from the multivariate dynamics of Eq. [4] under the 
approximation that stochastic variables (iV(r), N 2 (t), . . . , A^_i(r)) = (n, n 2 , 
. . . , Tif_i) are fixed in time such that = for i, j < /. In this scenario, 
the summations on the right-hand side of Eq. [4] may begin at / 

Given the identity 

-Pfi Pn-Pri2 \n ' ' ' Pnt—i |ii,ti2i-i n /-2 -^"/v! rt sl n ! n 2)---i rl /— 1 ' ( B- H) 
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the stationary factor P n P n2 \ n ■ ■ ■ i 3 n/_i|n,na,...,n/_a cancels from both sides of 
Eg. [BIO] to yield 



dP r 



f\n,n2,--,nf- 



dr 

s s 



^ ] ^ ] {Tijn+e i -ejPnf+ef i -efj\n,n2,...,n f _ 1 TjinPfif\n,n2,...,nf^ 1 ) ^iji 



i=f j=fJ¥=i 



0©< 

(B.12) 



where n/ = (n/, . . . , ns) and e/j is an (S — f + 1) -dimensional unit vector 
along the ith-direction. Now let 

Jl rif J^—n— 712 «/ JL—n—n'2 rij 

6= J2 ••• 2 5(J z -n-n 2 n s ), (B.13) 

n /+1 =0 "s=0 



so that 



nf\n,ri2,...,nf_ 



J L ny 



n 2 ,...,n / _ 1 - 



(B.14) 



Applying Eq. B.13 to both sides of Eq. |B.12 , we obtain 



dP r 



f\n,n2,~.,nf- 



dr 

- JL z* i i , 



PjinPnf\n,n2,...,nf^i^) ®ij 



Jl n f s 

+ E * E (w,-^«, + «-„-«-,>« 

i=/+i 



^j/^^%Kra 2 ,...,n/_i) @/j 



+ 4 f*t(W ta . 

i=/+l 



~PfinPfi[\ri,n2,...,nf—i) ®«/ - 

(B.15) 
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By inspection, the first term is identically zero and the remaining terms 
generate the right-hand side of Eq. B.7[ namely 

J Y'* i ^ W**. 

J=/+l 

^nj+l|n,ri2vi rl /— 1 ® ("^i (^/ ^))-^ri/+l|n.,n2v) n /-l ' 

J L re/ 5 

^ TjfnPrtf\n,n2,...,nf-x®fj 

9nf\n,n2,...,rif_ 1 

Q(J L - (n/ + l))Pn/|n,n 2 ,...,n/_ 1) 

T'^«*vw*. ^ 

i=f+l 

9nt— l\n,n2y,nf— i — I|re,n2, ...,n/_i j 

V' £ 

»=/+i 

7"n/|n,n2,.--,n/_i® l)-Fn/|n,ri2,---,n/_i - (B.16) 



The stationary distribution of Eq. B.7 is a Polya distribution (Johnson et al. 



1992) 



Dim* 
nf\n,ri2,---,n.f_ 



N, 



Pf 



J L - n,f-i\ B {(j) on Xf + n f ,J L - n/_i + on (1 - x/) - n/) 
n, j B 



l Xf, J L - n/_i + ort (1 - x/)) 



where 



(B.17) 



N 



B (Jl ~ Uf-l + 4> on (1 - Xf-i) ,(f) on (1 - X/)) 
Pf B (0 on (1 - X/_i) , J L - hf-l + <j)on (1 - X/))' 



(B.18) 



and 



^on lo 



l + n(w-l)/(J L -l) 
1 + x(w — 1) 



(B.19) 



Plugging Eqs. [21] and |B.17| into |B.6[ we obtain the approximate sampling 
distribution of Eq. [28} To validate Eq. [28l we demonstrate approximate 
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detailed balance, as denned by Eq. 10, for large- Jl nearly neutral commu- 
nities where S, w — 1 << Yld=i w i n i ~ 1 = Jl — 1 + n(w — 1) such that 
e = (w — 1)/(Jl — 1 + n{w — 1)) is a small number. For i,j > 2, detailed 
balance is exact. But for i — 1 and j > 2, we have 



pnn* 
n+ei— Sj 



n + 1 1 1 — m Q rij — 1 + 4> on Xj 
rij w 1 — m n + n x 

((1 - x)0 on+ i)j L _ w -i -pr (<j>onXk) nk ^ 5jk 

((1 



(B.20) 



where 



.l + n(«;-l)/(j£-l) 



1 + x(u» — 1) 
Given (a(l + e)) n ~ (a) n + O(e), we find 



pnn* 

n 

pnn* 



± ljn+e\—ej 



T - 



+ 0{Se). 



The case of i > 2 and j = 1 is similar. 



(B.21) 



(B.22) 



Appendix C. Asymptotic Expansions for Hypergeometric Func- 
tions 



Calculating Eqs. |2TJ [25j [27} [28j |A.8[ |A.11[ and |A. 14] for large communities 
requires computationally intensive evaluations of hypergeometric functions. 
To address this problem, one of us (N.M.T.) developed previously unknown 
asymptotic expansions. All required expansions are summarized here. Rel- 



Wong (2001); Gil et al. (2007) 



evant details can be found in Abramowitz and Stegun (1965); Luke (1969); 



Appendix C.l. Expanding 3-F 2 (l — Jl, 1 + A, 2; 2 — £, 1; rf) 
Using the reduction formula 



iF 2 (a,b,2;c, l;z) 



abz 



2 F 1 (a + 1, b + 1; c + 1; z) + 2 F 1 (a, b; c; z), (C.l.0.1) 



this case can be expanded with the methods of ^Appendix C.2 
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Appendix C.2. Expanding 2 -Pi (a — Jl, ol + A; a + 1 — 
Appendix C.2.1. Notation 
We write 

a = a — J L , b = a + (3 + {iJl, c = a + 7 + pJ^, (C.2.1.1) 

with a = 0, 1, 2 and J L a positive integer. In terms of w, m, x, and m we 
have 

Tfl T 

1 — m + a;(w — 1) 

and 

1 - xwm a + x(w - 1) 

7 = ^ —, tv- P = -7- (C.2.1. 3) 

The asymptotic behaviour will be considered of the Gauss hypergeometric 
function 

F = 2 F 1 (a,b;c;r ] ), (C.2.1.4) 

for large- Jl, where 

1 — m + x(w — 1) . _ n „ . 

77 = w- ^, C.2.1.5 

and 

we (0,oo), x,m,m G (0,1). (C.2.1. 6) 

Appendix C.2.2. The neutral case: w — 1, m — m 
In this case 

1 — mi „ „ _ . 

77 = 1, p=- , p=-- . C.2.2.1 

1 — m 1 — m 

The exact relation 

(c — b)„ r(c)T(c — b + n) ,^ 

2 F 1 (-n, b; c; 1 = = > K -(, n = 0,l,2,..., C.2.2.2 

(c) n T(c + n)T(c- b) 

can be used, together with the asymptotic estimate of the ratio of gamma 
functions 

E^±^ = n *-»(i + (i/ n )) > n ^oo. (C.2.2.3) 
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Appendix C.2.3. Critical values 

Considered as functions of w, \i and p become unbounded at w = w c ^ 
and w = w Cp , respectively, where 

m + x-1 l-x 
w c = , w Cp = . C.2.3.1 

x m Q — x 

The case w — >■ w c 

In this case r\ becomes small, b becomes unbounded, but the product brj 
remains finite. The kth term of the standard power series of F becomes (see 



also (C.2.2.3)) 



with 



where 



( a )k(b)k k ( Q )fc yk 
ti/ \ V ~ 777 \ z > (L.2.6.2) 
k\(c)k k!(c )k 



z= lim br] = u + vJi, cq = lim c = 70 + PoJl, (C.2.3.3) 



mx(m + x — 1) , _ 

u = 1 ; v = -u, C.2.3.4 

mx — m [m + x — 1) 

and 

x(m(l -m ) +m Q (l - x)) 

7o = 7 : 7-, — , Po = -7o- C.2.3.5 

mi — m {m + x — 1) 

It follows that F approaches a confluent hypergeometric function: 

2 F 1 (a,b;c;r]) ->■ ^(a; Co; (C.2.3. 6) 

Further action is needed to obtain an asymptotic approximation of the \F\- 
function. 

The case w — >■ w c „ 

In this case r\ and c become unbounded, but the ratio rj/c remains finite. 
The kth term of the standard power series of F becomes 

(a)k(b)k k ( Q )fcQo)fc , n o o 7 n 

ibi^ ~ (c - 2 - 3 - 7) 

with 

z = lim c/t] = u + vJl, b = lim b = (3 + p J L , (C.2.3. 8) 

w— >w Cp w— >w Cp 
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where 

m (m - x)(l - x) 

u = ■ — , v = -u, (C.2.3.9) 

mx — m [m + x — 1) 

and 

A> = ™*K>-s) p = -A, (C.2.3.10) 

ms — m [m + x — 1) 

It follows that .F approaches a 2-fo hypergeometric function 

a fi(a, b; c; ij) 2 F (a, b ; -; 1/*) = £ ^T^> (C.2.3.11) 

fc=0 

because a is a negative integer. This function can be expressed in terms of 
the Kummer ^/-function 

2 F (a, 6 ; -; l/z) = (-z) a U{a, 1 + a - 6 , -2). (C.2.3.12) 

Further action is needed to obtain an asymptotic approximation of the U- 
function. 

Appendix C.2-4- Expansion A 
An integral representation is 

2J F\(a, 6; c; r>) = T{b ^_ b) J ^(l ~ t) M (l ~ dt, (C.2.4.1) 

valid for c > b > 0, rj < 1. This integral can be used when p > p > 0, 77 < 1. 
As an example, consider 

r = 3, m=|, m„ = i x = |. (C. 2.4.2) 

This gives 

6 = a + ^(Jl - 1), c = a + 5(J L -l), p = jj, p = 5, 77 = -11. 

(C.2.4.3) 

In this case the integrand becomes small at t — and t = 1, and there is 
a maximum of the integrand at t — t\, with t\ G (0, 1). This point gives the 
main contribution. 

Write QC.2.4.1D as 

2 F 1 (a, b- c; 77) = ^77^ / t a+f3 -\l - ty~P-\l - t V )- a e- jL ^ dt, 



T(b)T(c-b) 



(C.2.4.4) 
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where 

<j){t) = -li ln(t) -{p-Li) ln(l -t)- ln(l - tr]). (C.2.4.5) 



The saddle points to and t\ are the zeros of 4>'{t). For the example (C.2.4.2) 
this gives 

t = -0.01169- •• , h = 0.1178- ■■ , (C. 2.4.6) 

and 

= -0.02136- •• , 0"(ti) = 35.83- •• . (C.2.4.7) 
An asymptotic approximation follows from the substitution 

<t>{t) - <j>(t x ) = ^"(ti)s 2 , sign(t - t x ) = sign(s), (C.2.4.8) 

which gives 

2 F x (a, 6; c; 7?) = —^L— e - J ^ f°° /( s ) e ~f ^"(^ 2 ds, (C.2.4.9) 

r(6)r(c-6) y_ oc 

where 

/(s) = r+^^l - ^t-^-^I - (C.2.4.10) 

Because locally at t = ti (or s = 0), t = ti + s + 0(s 2 ), we have dt/ds = 1 at 
s = 0, and 

/(0) = ^-'(l - tiy- p -\l - t x rj)- a . (C.2.4.11) 
This gives the first order approximation 

2 F 1 (a,b;c;r ] ) ~ r n3^_ b f ~ J ^f(0) f e^"^ ds, (C.2.4.12) 
that is 



(C.2.4.13) 
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Appendix C.2.5. Expansion B 

Another integral representation is 



2 Fi(a,b; c;r)) 



r(i + b - c) 
r(6)r(i - c) 



f-\t + l) c - b -\l + t7])- a dt, (C.2.5.1) 



which is only valid for a = 0, —1, —2, . . . and c < a + 1. It can be verified by 
expanding (1 + trj)~ a in powers of rj. 

We have p > and p < —1, and because (see (C.2.1.2), (C.2.1.3) and 

mx p 



(C.2.1.5)) 



V 



(1 — m x) p' 



(C.2.5.2) 



we see that r\ > 0. 

As an example, consider 



This gives 

b = a + 3(Jl — 1), c = a + - Jl), p = 3, p 



(C.2.5.3) 



.15 „ = JL 

13' ' 13' 

(C.2.5.4) 



Write (C.2.5.1) as 



2 F 1 (a,b; c;r)) 



r(i + 6- c ) 
r(Wi-c) 



1?+f*-\t + l) 7 "' 3 " 1 (l + tri)- a e- jLm dt, 

(C.2.5.5) 



where 



ifj(t) = -p ln(t) - (p - p) ln(t + 1) - ln(l + trj). 



The saddle points to and t\ are for the example (C.2.5.3) 



to = -74.89- •• , ti = 3.385 



and 



ipih) = 2.251 • • • , = 0.04951 

An asymptotic approximation follows from the substitution 

- = l^"(*i)s 2 , si ^{t ~ h) = sign(s), 



(C.2.5.6) 

(C.2.5.7) 
(C.2.5.8) 

(C.2.5.9) 
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which gives 

2Fi(a,b;c-,r ] ) = ^l + t ) ~ C \ e- J ^ f° g{s)e~\ J ^"^ ds, (C.2.5.10) 

r(6;r(i - c) 

where 

g {s) = t a+l3 -\i + ty-p-\i + tv)~ a< f- s (c.2.5.11) 

Because locally at t — t\ (or s — 0), t — t\ + s + 0(s 2 ), we have dt/ds = 1 at 
s = 0, and 

0(0) = % +p -\l + hy-^l + t lV )- a . (C.2.5.12) 
This gives the first order approximation 

2*i(a, b; c; V ) ~ ^j"* ~ c) e - ^> g(0) f° e -^"(^ 2 ^ (C .2.5.13) 

r(&)r(i-c) 

that is 



(C.2.5.14) 

Appendix C.2.6. Expansion C 

If fi < p < — 1 and ^ < 0, apply the transformation 

2 Fi(a, 6; c; rj) = (1 - r/^i^a, 6'; c; r/), (C.2.6.1) 

where 

b' = c-b = p , + t i , J L , = 1- P, fi' = p-fi, »7' = -^r- (C.2.6.2) 

77 — 1 

Now, 

// > 0, p < -1, 7/ > 0, (C. 2.6.3) 

and it follows that Expansion B, § Appendix C.2."5 , applies to the Gauss 
function on the right-hand side of (C.2.6.1). 
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Appendix C.2.7. General cases for all non-critical values 
1. w c , w c < 

For all w > 0, we have p > 0, p < — 1, and v > 0, so use Expansion 



B, ^Appendix C.2.5 



2. w c > 0, w c < 



For all to c > w > 0, we have p < —1, p < —1, and r\ < 0, so use 
Expansion C, § Appendix C.2.6 



For all w > w c , we have fx > 0, p < —1, and ^ > 0, so use 



Expansion B, ^Appendix C.2.5 



3. w c < 0, w c > 



For all w Cp > w > 0, we have p > 0, p < —1, and 7/ > 0, so use 



Expansion B, § Appendix C.2.5 



For all w > w Cp , we have p > p > and rj < 0, so use Expansion 
A, ^Appendix C.2.4 



4. uv„ > w r „ > 



For all w c > w > 0, we have /z < —1, p < —1, and i] < 0, so use 
Expansion C, § Appendix C.2.6 



For all w Cp > w > w c , we have p > 0, p < —1, and > 0, so use 
Expansion B, ^Appendix C.2.4 



For all w > w Cp , we have p > p > and r\ < 0, so use Expansion 
A, §Appendix C.2.4 



Appendix C.3. Expanding 2-F\(l — Jm, 1; 2 — ^a/; Vm) 
Appendix C.3.1. Notation 
We write 

a = 1 — J A ./, & = 1, c = cr + r J M , 

with 



a = l + 



1 



1 



(C.3.1.1) 



(C.3.1.2) 



1 — wv a 1 — WV Q 

The asymptotic behaviour will be considered of the Gauss hypergeometric 
function 

F = 2 F 1 (a,b;c;r ]M ) (C.3.1.3) 

for large- Jm, where 

(C.3.1.4) 



Vm 



and 



w e (o, oo), v, v Q e (o, l). 



(C.3.1.5) 
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Appendix C.3.2. The neutral case: w = 1, v = v 

In this case T)m — 1 an d (C.2.2.2) can be used to get an exact result in 
terms of gamma functions. 

Appendix C. 3. 3. The critical case w Cv =\jv 



In this case we have (see also £ Appendix C.2.3) 

2 F 1 {a,b;c;r ]M ) 2 F {a, b; -; l/z) = (-z) a U(a, 1 + a - b, -z), (C.3.3.1) 
where 

z = - (Jm ~ 1)z/ ° . (C.3.3.2) 
1 — v 

Appendix C.3.4- The case < w < w Cvo 
Use the integral representation 

^(a, b- c; VM ) = r ( 6)r(ilc) ^f®*''"*® dt > ( C - 3A1 ) 

where 

/(*) = (l+ty- b -\l + VM t)-\ <f>{t) = -Tln(l + t)-ln(l + VM t). (C.3.4.2) 
The saddle point to follows from solving = 0. This gives 

m = ~— — , to = ~ JL P^- (C3.4.3) 

l + t l + r] M t r)M{T + l) 

In terms of w and v 

(l-w)(l- TO + w) 
z/ (l — i/Jitr 

1. If < w < 1/(1 — v) = u> Ct ,, then the saddle point is negative, and we 
can substitute s = (pit), giving 

TYi + h _ r ) r°° 
2 F 1 (a,b;c; VM ) = -^- (/ s^g^e'^ ds, (C.3.4.5) 

r(6)r(i -c) j 

where 
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Apply Watson's lemma by expanding g(s) = Y^k=o9k sk ^° obtain 



2 -Fi(a,6; c; i] M ) 



r(l + &-c) y> r(6 + fc)ff fc 

r(6)r(i - c) ^ ' 



(C.3.4.7) 



To compute the coefficients we first expand t = ^2T=i ^kS k ■ The 
coefficients tk follow from inverting the expansion 

oo 

s = -r ln(l + t) - ln(l + 7] M t) = si = -r - 7? M . (C.3.4.8) 

This gives 



fe=i 



1 — wv 



T + T)M 1—W+WV 



(C.3.4.9) 



and for the first coefficient in the expansion (C.3.4.7) go = g(0) = t\. 
This gives 



2 Fi (a, b; c; r) M ) 



r(l + 6-c) / h 



V(l-c) 



J.. 



M 



(C.3.4.10) 



2. If w > w ClJ , then t is positive, and Laplace's method can be used, as 
in § §Appendix C.2.4~l |Appendix C.2.5| We substitute 



V 2 M (r + lf 



\<j>"{t )s 2 = <j>{t) - c/>(t ), 4>"(t ) - r( 1)2 ! 



(C.3.4.11) 



and obtain 



2 F 1 (a,b;c;r] M ) ~ 



2tt 



r(i + 6-c) 



r(6)r(i - c) y W(t ) 



-J M <t>(to)-ih-l 



(C.3.4.12) 



where t is given in (C.3.4.3). 



3. If w = w Cv , then to = and Laplace's method on a half-infinite interval 
can be used. 

Appendix C.3.5. The case w > w Cvo 
Use the integral representation 



2 F 1 (a,b; c; r) M ) 



Tie) 



T(b)T(c-b) 
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f{t)e~ jMm dt, 



(C.3.5.1) 



where 

/(t) = t b -\l - ty- b -\l - mi t)-\ <P(t) = -r ln(l - 1) - ln(l - VM t). 

(C.3.5.2) 

The saddle point to follows from solving <p'(t) = 0. This gives 

<p{t) = - - + - -, t = — - f — , 7T- C.3.5.3 

1 - t 1 - TJ M t 7] M (T + 1) 

In terms of w and v 

v [\ — v)w 2 

1. If w < w c „, then to < and Watson's lemma should be used. The 
result is 

2 F 1 (a,fo;c;7 7M )~-^-f|-') , tl = —±—. (C.3.5.5) 
T(c-b) \J M J Vm + t 

2. If w > w Cv , then the saddle point t is always inside the interval (0, 1), 
with t — > 1 if w — > oo. Laplace's method should be used. This gives 



T(c) 2vr 



(C.3.5.6) 

where /, and to are given in ( C.3.5.2[ )-( C.3.5.3) and 

3. If w = u> Ci/ then t = and Laplace's method on a half-infinite interval 
can be used. 

Appendix C.Jf.. Expanding 3-F 2 (l — Jm, 1,1;2, 2 — £m',Vm) 
Appendix CJ^.l. Notation 
We write 

a = l — J M) c = o- + tJ m , (C.4.1.1) 

with 

a = l + — 1 , r=-— 1 , (C.4.1.2) 

1 — wv a 1 — wv a 
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The asymptotic behaviour will be considered of the hypergeometric function 

(C.4.1.3) 



for large- J M , where 
and 

Note that 



F = 3 F 2 (a, 1, l;c,2;77 M ), 
w(l — u) 

Vm 



1 — wv 
re(0,oo), v, v Q G (0, 1). 



(C.4.1.4) 
(C.4.1.5) 



iF 2 (a, 1, 1; c, 2; = ^ J Vm 



k=0 



k\ (c) t (2), 



k=0 



(a)k Vm 
(c) k k + l 



(C.4.1.6) 



and because a is a negative integer this series terminates at k — —a. 



Appendix C.4-2. Representation in terms of a Laplace integral 
We substitute 



1 



jfe+1 



< k+l ^dw, 



(C.4.2.1) 



and obtain 



PCX) 

3 F 2 (a,l,l;c,2; %/ ) = / e -*' 2 F 1 (a,l;c,z)dw, z = nM e~ w . (C.4.2.2) 

Jo 

For this representation we use the results of § § Appendix C.3T4 Appendix 



C.3.5 


when, 


and ( 


C.3.5.3 



C.3.5 when, with tjm replaced by t\m& w , the saddle points to of (C.3.4.3) 



1. The case w < w Cu and w < w Cv . For this case we use the results in 



dC.3.4.5[ )-( [a3XT0| with 6 = 1. We have 

1 -c 



iFAa, l;c; z] 



J, 



M 



g + 0(1/ J M ) 



(C.4.2.3) 



where go = t\ and t\ = — l/(r + rjMe w ). This gives 



3 F 2 (a, 1,1;c,2;?7m) 



c- 1 



J 



M Jo 



dw 



t + Vm£~ 



+0(1/J M ). (C.4.2.4) 



Evaluating the integral we obtain 



3 F 2 (a, 1, 1; c, 2; tjm) = ° ±M± + Vm/t) + ^ (CA2 . 5 ) 



tJ, 



M 



Vm/t 
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2. The case w > w Cvo and w < w Cv . In this case we use the results in 
(C.3.5.1 )-(C.3.5.5), again, with 6 = 1 and t\m replaced with t]m£~ w ■ 
We use Watson's lemma for (C. 3. 4.11) by substituting s = <fi(t) and 
obtain 

:i - ty- 2 dt 

ds 



2 Fi(a,l;c; z) = (c - 1) 



-Jms 



g(s)ds, g(s) 



1 -te~ 

(C.4.2.6) 

Expanding (7 at s = we have g(s) = go+0(s), with g = l/(r+riMe~ w ) 
and, as in the above case, we derive 



iF 2 (a,l,l;c,2;r) M ) 



c - 1 ln(l + t] M /t) 



tJ m Vm/ t 
Appendix C.4-3. Summing the series by integration 



+ 0{1/J M ). (C.4.2.7) 



We replace the Pochhammer symbols in (C.4.1.6) by representations in 
terms of the gamma functions 



(«)* 



:-r. 



■ r(i-a) 

r(l — a — k)' 



(C.4.3.1) 



and replace the gamma functions with large positive argument by their 
asymptotic forms that follow from 



- az (az) az+b -k 



z — Y 00, a > 0. 



(C.4.3.2) 



r(az + 6) ~ v 27r e 
This gives the remaining two cases. 

1. T/ie case w < w Cvo and w > w Cu . In this case r < — 1 and t\m > 0. We 
replace the Pochhammer symbols in (C.4.1.6) with the second form in 

(C.4.3.3) 



QC.4.3.1). This gives 



where 



F(k) 



3 F 2 (a,l,l;c,2;r ?M ) = ^F(fc), 

k=0 

r(l-a) r, k M r(l-c-fc) 
r(l-c)Jfe + ir(l-a-A;)' 



(C.4.3.4) 



and we replace the summation in (C.4.3.3) by integration, invoking 
Euler's summation formula, or the compound trapezoidal rule 



3 F 2 (a, 1, 1;c,2;t]m) 



T(l-a) 

fa - c) 



■hi 



v k M r(i-c-fc) 

k + l T(l - a- k) 



dk. (C.4.3.5) 
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In Euler's summation formula additional terms occur but in the present 
case they can be neglected. 

We replace the gamma functions by their asymptotic estimates follow- 



ing from (C. 4.3.2). This gives 



. T(l - a) [ Jm y/1 — a — k „ ln 

3 F 2 (a, 1, 1; c, 2; m ) ~ ( / V = dfc, C.4.3.6 

r(i - cj Jo fe + l vi - c - k 

where 

0(fc) = — khiriM — (1 — c — k) ln(l — c — /c) + (1 — a — k) ln(l — a — k). 

(C.4.3.7) 

Then, 

<j/(k) = -ln?7 M + hi(l-c- k) -ln(l - a- fc)), (C.4.3.8) 
and the derivative vanishes for k — k s , where 



c-l + rj M {l-a) a -I 

k s = ; = k + k\JMi k 



1 - Vm 



Vm - 1 



Vm -I' 
(C.4.3.9) 



and < &i < 1. The dominant point of the integral in (C.4.3.6) is 
k = k s , and we apply Laplace's method. We substitute 

U"{k s )s 2 = m - <p(k s ), 4>"(k s ) = - + J2r). 



JmVm(t + 1) 



(C.4.3.10) 



This gives 



jF 2 (l - J M , 1, 1; a + r Jfcf, 2; r/ M ) 



2tt 



F(k s ), (C.4.3.11) 



where -F(fc) is given in (C.4.3.4). After using (C.4.3.2) we obtain 



jF 2 (l - J M , 1, 1; a + tJ m , 2; r/ M ) 



2tt - 1 



<P"{k s ) (rj M + t)J m \ ^~ ' V 1 - 



l + r 



(C.4.3.12) 
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2. The case w > w Cvg and w > w Cv . In this case r > and rj M < 0. We 
use (C.4.1.6), replacing the Pochhammer symbol (a)k by the second 



form of (C.4.3.1) and (c)& by the first. This gives 



where 



F(k) 



F 2 (a,l,l;c,2;ri M ) = J2 F (k): 
T(J M )T(c)(- VM ) k 



(C.4.3.13) 



(C.4.3.14) 



(k + l)T(J M -k)T(c + ky 
and we replace the summation in ( C.4.3.13"] ) by integration, 

3 F 2 (a,lA;c,2;r]M) ~ / F{k)dk. (C.4.3.15) 

J o 



Applying the asymptotic estimates of gamma functions in (C.4.3.2), we 
obtain 

m = T{J 2n(k C +^ JM V(JM ~ k)(C + k)e ~ m > (CA3 " 16) 

where 

4>{k) = (c + k)ln(c + k) + (J M -k)\n(J M -k)-kln(-T] M ). (C.4.3.17) 
We have 

4>'(k) = ln(c + k) - ln(J M - k) - H- VM ), (C.4.3.18) 
and the saddle point given by 



C + T] M J M 

k s = — = k + kiJ M , k 



cr , r + riM 



Vm - 1 
(C.4.3.19) 

where, again, < ki < 1. The dominant point of the integral in 
(C.4.3.15) is k = k s , and we apply Laplace's method to this integral. 
We substitute 

l<p"(k s )s 2 = <f>(k) - <p(k s ), <f>"(k s ) = - r^-r^\ + (V^m), 

JmVm{t + 1) 

(C.4.3.20) 
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where 4>{k) is given in (C.4.3.17). This gives 
3 F 2 (1 - J M , 1, 1; a + tJ m , 2; T] M ) r - 



2tt 



F{k s ), (C.4.3.21) 



where F(k) is given in (C.4.3.16). After using QC.4.3.2D we obtain, as 
in the case above, 

3 F 2 (1 - J M , 1, 1; a + rJ M , 2; rj M ) 



2?r t)m - 1 (_Vm_\ ^ c ( J: 



4>"(k s ) (rjM + r)J M V t ) \l-r} M 



(C.4.3.22) 
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